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Abstract 

We carry out a systematic study of entanglement entropy in relativistic quantum field theory. 
This is defined as the von Neumann entropy Sa = — Tr PA log PA corresponding to the reduced 
density matrix pa of a subsystem A. For the case of a 1+1-dimensional critical system, whose 
continuum limit is a conformal field theory with central charge c, we re-derive the result Sa ~ 
(c/3) log t of Holzhey et al. when A is a finite interval of length I in an infinite system, and extend 
it to many other cases: finite systems, finite temperatures, and when A consists of an arbitrary 
number of disjoint intervals (See note added). For such a system away from its critical point, 
when the correlation length ^ is large but finite, we show that Sa ~ .A(c/6)log£, where A is the 
number of boundary points of A. These results are verified for a free massive field theory, which is 
also used to confirm a scaling ansatz for the case of finite-size off-critical systems, and for integrable 
lattice models, such as the Ising and XXZ models, which are solvable by corner transfer matrix 
methods. Finally the free-field results are extended to higher dimensions, and used to motivate a 
scaling form for the singular part of the entanglement entropy near a quantum phase transition. 
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I. INTRODUCTION. 



Recently there has been considerable interest in formulating measures of quantum entan- 
glement and applying them to extended quantum systems with many degrees of freedom, 
such as quantum spin chains. 

One of these measures |l| is entanglement entropy. Suppose the whole system is in a pure 
quantum state with density matrix p = |^ r )(^ f |, and an observer A measures only a 
subset A of a complete set of commuting observables, while another observer B may measure 
the remainder. A's reduced density matrix is pa = Tr# p- The entanglement entropy is just 
the von Neumann entropy Sa = — Tr^ pa log pa associated with this reduced density matrix. 
It is easy to see that Sa = Sb- For an unentangled product state, Sa = 0. Conversely, Sa 
should be a maximum for a maximally entangled state. 

For example, for a system with two binary (spin--) degrees of freedom, with \^f) = 
cos 6*| tl) + sin#| where A observes only the first spin and B the second, Sa takes its 
maximum value of log 2 when cos 2 # = |, which agrees with our intuitive idea of maximal 
entanglement. For this system it has been shown, [1] even in the partially entangled case, 
that if there are M copies of the state available, by making only local operations A can 
produce M' < M states which are maximally entangled. The optimal conversion ratio 
M'/M is given, for large M, by S A -[28\ 

Although there are other measures of entanglement J2| the entropy is most readily suited 
to analytic investigation. In several papers jj, 0, Q], the concept has been applied to 
quantum spin chains. Typically, the subset A consists of a commuting set of components 
of the spin degrees of freedom in some interval of length £, in an infinitely long chain. It is 
found that the entanglement entropy generally tends to a finite value as I increases, but that 
this value diverges as the system approaches a quantum critical point. At such a critical 
point, the entropy grows proportional to log^ for large I. 

Close to a quantum critical point, where the correlation length £ is much larger than the 
lattice spacing a, the low-lying excitations and the long-distance behaviour of the correlations 
in the ground state of a quantum spin chain are believed to be described by a quantum field 
theory in 1+1 dimensions. If the dispersion relation of the low-lying excitations is linear 
for wave numbers k such that <C \k\ <C a -1 , the field theory is relativistic. We shall 
consider only those cases in this paper. At the critical point, where £ _1 = 0, the field theory 
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is massless, and is a conformal field theory (CFT). 

In this case, the von Neumann entropy of subsystem corresponding to an interval of 
length £ was calculated some time ago by Holzhey et al.[z], in the context of black hole 
physics (although that connection has been questioned), where it was termed 'geometric' 
entropy. Using methods of conformal field theory, based in part on earlier work of Cardy 
and Pescheljs], they found Sa ~ (c/3) \og(£/a), where c is the conformal anomaly number 
(sometimes called the central charge) of the corresponding CFT. 

This result has been verified by analytic and numerical calculations on integrable quantum 
spin chains corresponding to CFTs with c = | and c = 1. j^l, Q, Q] 

In this paper, we first put the CFT arguments of Holzhey et al. [?| on a more systematic 
basis, and generalise their result in a number of ways. Our methods are based on a formula 
for the entropy in terms of the partition function in the path integral formulation of the 
quantum theory as a euclidean field theory on an n-sheeted Riemann surface, in the limit 
n — > 1. 

For a 1+1-dimensional theory at a critical point, we derive analogous formulas for the 
entropy in the cases when the subsystem A consists of an arbitrary number of disjoint 
intervals of the real line (see Eq. f¥Tj) . and when the whole system has itself a finite length 
L. For example, for the case when A is a single interval of length £, and periodic boundary 
conditions are imposed on the whole system, we find 

S A = (c/3) log ((L/na) sm(n£/L)) + c[ . (1) 

On the other hand, for a finite system of total length L with open boundaries, divided at 
some interior point into an interval A of length £ and its complement, we find 

Sa = (c/6) log ((L/vra) sm(n£/L)) + 2g + c[ , (2) 

n 

where g is the boundary entropy of Affleck and Ludwig[9j. We also treat the case when the 
system is infinitely long but is in a thermal mixed state at finite temperature f3~ l : 

S A = (c/3) log (GS/tto) smh(n£/(3)) + c[ . (3) 

In all these cases, the constant c[ is the same, but non-universal. 

For a massive 1+1-dimensional relativistic QFT (which corresponds to an off-critical 
quantum spin chain where the correlation length £ ^> a) the simplest results are for an 



infinite system divided at some point into two semi-infinite pieces. In this case we verify 
that the entanglement entropy is finite, and derive the universal formula 

S A ~ (c/6) logtf/a) . (4) 

In the more general case when A consists of a collection of disjoint intervals, each of length 
^> £, we expect (@| to be multiplied by a factor A which counts the number of boundary 
points between A and B (the Id analogue of surface area.) 

For A — 1 the entropy is exactly calculable in the case of a free field theory. We verify 
the above formula, and exhibit the finite-size cross-over which occurs when £ is of the size 
L of the system. For a lattice model in this geometry, with L infinite, we point out that pa 
is simply related to Baxter's corner transfer matrix, and thus, for integrable models whose 
weights satisfy a Yang-Baxter relation, all its eigenvalues can be determined exactly. We 
treat explicitly the case of the Ising model and its anisotropic limit, the transverse Ising 
spin chain, and also the XXZ model, computing exactly at all values of the coupling the 
finite part of the entropy S. This agrees with our continuum result (j3J) for c = |, 1 when 
the correlation length £ is large. 

The analysis for the free theory is straightforward to extend to higher dimensions, at least 
in suitable geometries. This leads to the well-known law, first found by Srednicki[10|], that 
the entropy should be proportional to the surface area A of the subsystem A. As pointed 
out by Srednicki, the coefficient of this term, for d > 1, depends on the UV cut-off and is 
therefore non- universal. However, based on our calculations, we propose that there should 
be a non-leading piece in Sa/A, proportional to £~( c( ^ 1 ) ) which depends in a singular way 
on the couplings near a quantum phase transition, and whose form is, moreover, universal. 

The layout of this paper is as follows. In the next section, we discuss the entropy in terms 
of the euclidean path integral on an n-sheeted Riemann surface. In Sec. IIHI we consider 
the 1+1-dimensional conformal case. We use the powerful methods of CFT to show that 
the partition function on the Riemann surface of interest is given, up to a constant, by a 
calculable correlation function of vertex operators in a c = 1 CFT in the complex plane. 
Similar results apply to a system with boundaries. This general result allows us to derive 
all the special cases described above. In Sec. [TV] we consider the case of a massive 1 + 1- 
dimensional field theory. We derive the result S ~ (c/6) log(£/a) from completely general 
properties of the stress tensor in the relevant geometry. This is then verified for a free bosonic 
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massive field. In the last part of this section we relate the lattice version of this problem to 
the corner transfer matrix, and compute the off-critical entropy for the case of the Ising and 
XXZ Heisenberg spin chains. In Sec. |V]we study the off-critical case in a finite-size system, 
propose a scaling law, and verify it for the case of a free massive field theory. We compute 
the relevant scaling function in a systematic sequence of approximations. Sec. I VII is devoted 
to the discussion of higher dimensions. 

While all this work was being carried out, some other related papers have appeared 
in the literature. In Ref. the result Sa ~ ( c /3) log(£/a), first found by Holzhey et 



al.jz], was obtained by the following argument: it was assumed that the entropy should be 
conformally invariant and therefore some function F(x) of the variable x — ((3 /it) sinh(7r£//3); 
by comparing with known case £ 3> (3 [12I . \\\ it was observed that F(x) should behave as 
(c/3) log x as x — > 00; finally it was assumed (with no justification being given) that the 
log x form for F is valid for all values of x; finally the limit (3 3> I was taken. In the present 
paper, we should stress, we have derived all these statements from first principles of CFT. 

Very recently Casini and Huertajy] have considered the case of two intervals A and B, 
and argued that the quantity F(A, B) = S(A) + S(B) - S(A n B) - S(A U B) is UV 
finite as a — > 0, and is given by a universal logarithmic function of the cross-ratio of the 
four end points. This corresponds to, and agrees with, our case N = 2. These authors, 
however, assume the conformal invariance of the entropy, while, once again, we stress that 
in the present paper we derive this from fundamental properties of the stress tensor. These 
authors also give a very nice alternative derivation of Zamolodchikov's c-theoremfrJ] based 
on this quantity F(A, B). 



II. VON NEUMANN ENTROPY AND RIEMANN SURFACES. 

Consider a lattice quantum theory in one space and one time dimension, initially on the 
infinite line. The lattice spacing is a, and the lattice sites are labelled by a discrete variable 
x. The domain of x can be finite, i.e. some interval of length L, semi-infinite, or infinite. 
Time is considered to be continuous. A complete set of local commuting observables will 
be denoted by {<p(x)}, and their eigenvalues and corresponding eigenstates by {4>{x)} and 
® x \{4>{x)}) respectively. For a bosonic lattice field theory, these will be the fundamental 
bosonic fields of the theory; for a spin model some particular component of the local spin. 
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The dynamics of the theory is described by a time-evolution operator H. The density matrix 
p in a thermal state at inverse temperature (3 is 

piWrmix')'}) = Z(P)-\{<f>(xy}\e^ 6 \{4>(x'y}) , (5) 

where Z(f3) = Tie~^ H is the partition function. 

This may be expressed in the standard way as a (euclidean) path integral: 

p = Z- 1 [ [d<t>{x, t)\ 1] 5(<j>{x, 0) - II K<t>{x, /3) - 0(x")") e~ s * , (6) 

J X X 

where Se = Jo LeAt, with Le the euclidean lagrangian. (For a spin model this would be 
replaced by a coherent state path integral.) 

The normalisation factor of the partition function ensures that Tr p = 1, and is found by 
setting {(f)(x)"} = {<p(x)'} and integrating over these variables. This has the effect of sewing 
together the edges along r = and r = (3 to form a cylinder of circumference 3. 



Now let A be a subsystem consisting of the points x in the disjoint 29| intervals 
(ui,vi), . . . , (un,vjst). An expression for the the reduced density matrix pa may be found 
from ([6]) by sewing together only those points x which are not in A. This will have the effect 
of leaving open cuts, one for each interval (uj,Vj), along the the line r = 0. 

We may then compute Tr p n Al for any positive integer n, by making n copies of the above, 
labelled by an integer k with 1 < k < n, and sewing them together cyclically along the the 
cuts so that 4>{x)' k = 4>(x)'l +1 (and </>(x)' n = <j>(x)") for all x G A. Let us denote the path 
integral on this n-sheeted structure by Z n (A). Then 

Trp" A = ^. (7) 

Now, since Trp^ = J2x^ n , where A are the eigenvalues of pa (which lie in [0, 1),) and 
since Tr pa = 1, it follows that the left hand side is absolutely convergent and therefore 
analytic for all Ren > 1. The derivative wrt n therefore also exists and is analytic in the 
region. Moreover, if the entropy Pa = — J2x A log A is finite, the limit as n —>■ 1+ of the first 
derivative converges to this value. 

We conclude that the right hand side of (j7j) has a unique analytic continuation to Re n > 1 
and that its first derivative at n = 1 gives the required entropy: 

n-+i dn Z n v ' 
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(Note that even before taking this limit, (JTj) gives an expression for the Tsallis entropy 15] 
(Trpl-l)/(l-n).) 

So far, everything has been in the discrete space domain. We now discuss the continuum 
limit, in which a — > keeping all other lengths fixed. The points x then assume real values, 
and the path integral is over fields 4>{x, r) on an n-sheeted Riemann surface, with branch 
points at Uj and Vj. In this limit, Se is supposed to go over into the euclidean action for a 
quantum field theory. We shall restrict attention to the case when this is Lorentz invariant, 
since the full power of relativistic field theory can then be brought to bear. The behaviour 
of partition functions in this limit has been well studied. In two dimensions, the logarithm 
of a general partition function Z in a domain with total area A and with boundaries of total 
length £ behaves as 

log Z = f ± Aa' 2 + f 2 Ca- 1 + ... (9) 

where f\ and fi are the non-universal bulk and boundary free energies. Note, however, that 
these leading terms cancel in the ratio of partition functions in (J7J). However, as was argued 
by Cardy and PeschelQ] there are also universal terms proportional to log a. These arise 
from points of non-zero curvature of the manifold and its boundary. In our case, these are 
conical singularities at the branch points. In fact, as we shall show, it is precisely these 
logarithmic terms which give rise to the non-trivial dependence of the final result for the 
entropy on the short-distance cut-off a. For the moment let us simply remark that, in order 
to achieve a finite limit as a — > 0, the right hand side of (0) should be multiplied by some 
renormalisation constant Z(A,n). Its dependence on a will emerge from the later analysis. 



III. ENTANGLEMENT ENTROPY IN 2D CONFORMAL FIELD THEORY. 

Now specialise the discussion of the previous section to the case when the field theory 
is relativistic and massless, i.e. a conformal field theory (CFT), with central charge c, and 
initially consider the case of zero temperature. 

We show that in this case the ratio of partition functions in ([7j) is the same as the corre- 
lation function arising from the insertion of primary scaling operators <& n (iij) and <$>_ n (vj), 
with scaling dimensions X n = 2A n = (c/12)(l — 1/n 2 ), into each of the n (disconnected) 
sheets. Moreover, this 2iV-point correlation function is computable from the Ward identities 
of CFT. 
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In the language of string theory, the objects we consider are correlators of orbifold points 
in theories whose target space consists of n copies of the given CFT. We expect that some 
of our results may therefore have appeared in the literature of the subject. 



A. Single interval 

We first consider the case N — 1 and no boundaries, that is the case considered by 
Holzhey et al. [3] of a single interval of length t in an infinitely long Id quantum system, at 
zero temperature. The conformal mapping w — > ( = (w — u)/(w — v) maps the branch points 
to (0, oo). This is then uniformised by the mapping ( — > z = £ 1//n = ((w — u)/(w — v)) 1 ^ 1 . 
This maps the whole of the n-sheeted Riemann surface lZ n to the z-plane C. Now consider 
the holomorphic component of the stress tensor T(w). This is related to the transformed 



stress tensor T(z) by[16j] 



T(w) = (dz/dwf T{z) + y|{z, w} , (10) 

where {z,w} is the Schwartzian derivative (z"'z' — |z" 2 )/z' 2 . In particular, taking the 
expectation value of this, and using (T(z))c = by translational and rotational invariance, 
we find 

irpf U C, , C(l-(l/n) 2 ) {V-Uf 

(T{w))n n = ^{z, w} = — — - . 11 

12 24, (w — u) z (w — v) z 



This is to be compared with the standard form[16|] of the correlator of T with two primary 
operators $ n (ii) and $_ n (u) which have the same complex scaling dimensions A n = A n = 
(c/24)(l-(l/n) 2 ): 

(T(w)$ n (u)^ n (v)) c = An s , (12) 

(w — u) z {w — V) Z [V — u) Z£An z (v — U) z ^ n 

where $± n are normalised so that ($ n (ti)$_ n (w))c — \v — M |- 2A «- 2A ™. ( [121) is equivalent to 
the conformal Ward identity [l(]] 

(T( W )<S> n (u)^ n (v)) C = (-r-^yi + + + —4-) (M^-n(v))c ■ 

\{w — u) z (w — v ) z w — uou w — V ov J 

(13) 

In writing the above, we are assuming that w is a complex coordinate on a single sheet C, 

which is now decoupled from the others. We have therefore shown that 

_ S[d4]T{w)e- s *<™ _ (T(w)<5> n (u)^ n (v)) c 
\ W/w.- m]e s E (n n ) ($ n ( w )d>_ n („)) c • 1 } 
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Now consider the effect of an infinitesimal conformal transformation w — > w' = w + a(w) of 
C which act identically on all the sheets of TZ n . The effect of this is to insert a factor 

/ a(w)T(w)dw / a(w)T(w)dw (15) 

2ni Jc 2m Jc 

into the path integral, where the contour C encircles the points u and v. The insertion of 
T(w) is given by (fl2"l) . Since this is to be inserted on each sheet, the right hand side gets 
multiplied by a factor n. 

Since the Ward identity (fT3l) determines all the properties under conformal transforma- 
tions, we conclude that the renormalised Z n (A)/Z n oc Trp^ behaves (apart from a possible 
overall constant) under scale and conformal transformations identically to the nth power 
of two-point function of a primary operator $ n with A n = A n = (c/24)(l — (1/n) 2 ). In 
particular, this means that 

Trp n A = c n ((v-u)/a)-^ n - 1 / n \ (16) 

where the exponent is just AnA n . The power of a (corresponding to the renormalisation 
constant Z) has been inserted so as the make the final result dimensionless, as it should 
be. The constants c n are not determined by this method. However c\ must be unity. 
Differentiating with respect to n and setting n = 1, we recover the result of Holzhey et al. 

The fact that Tr p n A transforms under a general conformal transformation as a 2-point 
function of primary operators $± n means that it is simple to compute in other geometries, 
obtained by a conformal mapping z — > z' = w(z), using the formula [3] 

Z2) . . .) = J] \w\zj) | 2A " ($(«>!, Wi)$(v)2, W 2 ) . . .) . (17) 

3 

For example, the transformation w — > w' = (/3/2tt) logw maps each sheet in the w-plane 
into an infinitely long cylinder of circumference j3. The sheets are now sewn together along 
a branch cut joining the images of the points u and v. By arranging this to lie parallel to 
the axis of the cylinder, we get an expression for Tr p\ in a thermal mixed state at finite 
temperature This leads to the result for the entropy 

S A (J3) ~ (c/3) log (((3/na) sinh(vr£//?)) + c[ . (18) 

For £ (3 we find Sa ~ ( c /3) \og{£/a) as before, while, in the opposite limit I ^> (3, 
Sa ~ ( 7r c/3)(£//3). In this limit, the von Neumann entropy is extensive, and its density 



agrees with that of the Gibbs entropy of an isolated system of length £, as obtained from 



the standard CFT expression 
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13| f3F ~ —(ttc/6)(£/0) for its free energy. 



Alternatively, we may orient the branch cut perpendicular to the axis of the cylinder, 
which, with the replacement (3 — > L, corresponds to the entropy of a subsystem of length 
£ in a finite Id system of length L, with periodic boundary conditions, in its ground state. 
This gives 

S A ~ (c/3) log ((L/na) sm(rt£/L)) + c[ . (19) 
Note that this expression is symmetric under i —>■ L — i. It is maximal when £ = L/2. 



B. Finite system with a boundary. 

Next consider the case when the Id system is a semi-infinite line, say [0, oo), and the 
subsystem A is the finite interval [0, £). The n-sheeted Riemann surface then consists of 
n copies of the half-plane x > 0, sewn together along < x < £, r = 0. Once again, we 
work initially at zero temperature. It is convenient to use the complex variable w = r + ix. 
The uniformising transformation is now z = {{w - i£)/{w + i£)) 1/n , which maps the whole 
Riemann surface to the unit disc \z\ < 1. In this geometry, (T(z)) = by rotational 



invariance, so that, using (flOl) . we find 

(TH)n n = ~ ( ^ , (20) 

(w — i£y\w + %£y 

where A n = (c/24)(l — n~ 2 ) as before. Note that in the half-plane, T and T are related by 
analytic continuation: T(w) = T(w)*. fl20l) has the same form as (T(w)Q n (i£)) , which follows 



from the Ward identities of boundary CFT 17]], with the normalisation ($„(?£)) = (2£) An . 



The analysis then proceeds in analogy with the previous case. We find 

Trp^ ~ 5 n {2£/a) (c/12)(n - 1/n) , (21) 

so that S A ~ (c/6) log(2£/a) + c\. 

Once again, this result can be conformally transformed into a number of other cases. At 
finite temperature we find 

S A (P) ~ (c/6) log ((/3/tto) sinh(2vr£//?)) + c[ . (22) 

By taking the limit when £ S> (3 we find the same extensive entropy as before. However, we 
can now identify c[ — c[ as the boundary entropy g, first discussed by Affleck and LudwigP 

10 



For a completely finite Id system, of length L, at zero temperature, divided into two 
pieces of lengths £ and L — £, we similarly find 



S A = (c/6) log ((L/tto) sin(vr£/L)) + 2g + c[ . 



(23) 



C. General case. 

For general N, the algebra is more complicated, but the method is the same. The 
uniformising transformation now has the form z = Y\i(w — Wi) a % with J2i a % — (so there is 
no singularity at infinity.) Here Wi can be Uj, Vj, or ±iuj in the case of a boundary. In our 
case, we have a, = ±l/n, but it is interesting to consider the more general transformation, 
and the notation is simpler. Once again we have 



(24) 



Consider {z,w}/z. As a function of w, this is meromorphic, has a double pole at each 
w = Wi, and is 0{w~ 2 ) as w — > oo. Hence it has the form 



(w - Wif 



Bi 



W — Wi 



(25) 



where J2iBi = 0. In order to determine Ai and B iy we need to compute z', etc, to second 
order in their singularities at Wi. After some algebra, we find 



+ E 



ft; 



w — Wi ^r. w — Wj 



{W - Wi) 2 W-Wi^Wi- Wj 



2oti — 3a 2 + af —3a>i + 3a 2 aj 

(W - Wi) 3 (W - Wi) 2 Wi - Wj 



+ 



Let 



a i 



■ ,■ Wi — Wj 



(26) 
(27) 
(28) 

(29) 



Then the coefficient of (w — Wi) is 



{[a;(a:j - l)(«i - 2) + 3aj(aj - l)Cj(w — w t ) H ] [a, + C^w - w,j) + 

-§ [a* (a* - 1) + 2a;C;(u> — Wi) H ] 2 } / [a, t + Cj(u> — iUj) H ] 2 , 
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from which we find after a little more algebra that 



I -a 2 
r>i = Oj . 

on 



Thus we have shown that 



(TH) 



12 



E 



|(i -^ 2 ) 

(w - Wi) 2 



4 i^^E 



a,- 



. , . Wi — Wj ] W — Wi 



This is to be compared with the conformal Ward identity 



(twii«=e 



+ 



d 



(w — Wi) 2 ' w — Wi dwi 
For these to be equivalent, we must have Aj = |(1 — a 2 ) and 



(II 

k 



E 



a,- 
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A necessary and sufficient condition for this is that 



io g <n**K)>- 

k 



(30) 
(31) 

(32) 
(33) 
(34) 



d 



a. 



dwi 



on 



E 



a 4 



Wi — Wj 



d 

dwi 



a k \ ~> a j 

Oik W k - Wj 



for each pair (i, k). This reduces to 



l-a 2 



Oik 



\-a\ 



OLi 



(35) 



(36) 



di (w k - Wi) 2 a k (wi - w k ) 2 ' 

that is, ctj = ±a k for each pair (i, k). Since J2i a i — the only way to satisfy this is to have 
ctj = affj, with <7j = ±1, and half the <7j = +1 and the remainder —1. Interestingly enough, 
this is precisely the case we need, with a — 1/n. 
If these conditions are satisfied, 



_d_ 

dw, 



iog(n^K)) = ^(i-« 2 )E 



ViOk 



r. Wi - w k 



(37) 



3 ' "' k^i 

so that 

log(n^K)} = ^(l-a 2 )^2^klog{wi-w k ) + E, (38) 

j k<i 

where E is independent of all the Wi. In the case with no boundary, E can depend however 
on the Wi. A similar calculation with T then gives a similar dependence. We conclude that 
Tr p n A behaves under conformal transformations in the same way as 



<n $ ,> <* iik - w^-^^iwi - w k )^~ a2 ^ 

3 3<k 



(39) 
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Taking now a = 1/n, and a = ±1 according as Wi = Uj or Vj, we find 

T n „ c * ( n J<fc K-^)(^-^) \ (c/6)(n " 1/n) (40) 

The overall constant is fixed in terms of the previously defined c n by taking the intervals to 
be far apart from each other, in comparison to their lengths. 

Differentiating with respect to n and setting n — 1, we find our main result of this section 

S A = 3 I lo S (Ofc - u j)/ a ) ~ lo § d u k ~ u j)/ a ) ~ lo S (K ~ v i)l a ) + Nc 'i ■ ( 41 ) 

\j<k j<k j<k J 

A similar expression holds in the case of a boundary, with half of the u>, corresponding 
to the image points. 

Finally we comment on the recent result of Casini and Huerta[G], which corresponds to 
N = 2. In fact, it may be generalised to the ratio of Tsallis entropies: from (HUj) we find 

S n (A)S n (B) = ( {v 1 -u 1 ){v 2 ~u 2 ) \ (c/6)( - 1/n) 
S n (AnB)S n (AuB) U« 2 -Wi)(^-^i)y ' 1 ] 

where S n (A) = Tr p\, and the expression in parentheses is the cross-ratio 7] of the four points. 
Notice that in this expression the dependence on the ultraviolet cut-off a disappeared, so 
have the non-universal numbers c n . Differentiating with respect to n gives the result of 
Casini and HuertaQ], who however assumed that the result should depend only on rj. 



IV. ENTROPY IN NON-CRITICAL 1+1-DIMENSIONAL MODELS 
A. Massive field theory - general case 

In this section we consider an infinite non-critical model in 1+1-dimensions, in the scaling 
limit where the lattice spacing a — > with the correlation length (inverse mass) fixed. This 
corresponds to a massive relativistic QFT. We first consider the case when the subset A 
is the negative real axis, so that the appropriate Riemann surface has branch points of 
order n at and infinity. However, for the non-critical case, the branch point at infinity is 
unimportant: we should arrive at the same expression by considering a finite system whose 
length L is much greater than £. 

Our argument parallels that of Zamolodchikov[14] for the proof of his famous c-theorem. 
Let us consider the expectation value of the stress tensor T^ u of a massive euclidean QFT 
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on such a Riemann surface. In complex coordinates, there are three non-zero components: 
T = T zz , T = Tzzi and the trace 9 = 4T zS = 4Tg z . These are related by the conservation 
equations 

d g T+\d z Q = 0; (43) 
d s T+\d s e = 0. (44) 

Consider the expectation values of these components. In the single-sheeted geometry, 
(T) and (T) both vanish, but (0) is constant and non-vanishing: it measures the explicit 
breaking of scale invariance in the non-critical system. In the n-sheeted geometry, however, 
they all acquire a non-trivial spatial dependence. By rotational invariance about the origin, 
they have the form 

(T{z,z)) = F n (zz)/z 2 ; (45) 
(9(3,;*)) -(6)! = G n {zz)/(zz); (46) 
(T(z,z)) = F n {zz)/z 2 . (47) 

From the conservation conditions (143!) we have 

(zz) (F' n + \G' n ) = \G n . (48) 

Now we expect that F n and G n both approach zero exponentially fast for \z\ ^> £, while 
in the opposite limit, on distance scales <C ^, they approach the CFT values (see previous 
section) F n -> (c/24)(l - n~ 2 ), G n -> 0. 

This means that if we define an effective C-function 

C n (R 2 ) = (F(R 2 ) + ±G(R 2 )) , (49) 

then 

R 2 Q^T ) Cn(R 2 ) = \G n (R 2 ). (50) 

If we were able to argue that G n < 0, that is (9)„ < (O)i, we would have found 
alternative formulation of the c-theorem. However, we can still derive an integrated form of 
the c-theorem, using the boundary conditions:}^] 

^^^ 2 ) = -(c/6)(l-^ 2 ), (51) 
14 



or equivalently 

J ((6> n - (6)0 d 2 R = -7m(c/6)(l - n- 2 ) , (52) 

where the integral is over the whole of the the n-sheeted surface. Now this integral (multi- 
plied by a factor 1/2-7T corresponding to the conventional normalisation of the stress tensor) 
measures the response of the free energy — log Z to a scale transformation, i.e. to a change 
in the mass m, since this is the only dimensionful parameter of the renormalised theory. 
Thus the left hand side is equal to 

- (2tt) m(d/dm) [log Z n -n log Z) , (53) 

giving finally 

= a™) (c/12)K1/n) , (54) 

where c n is a constant (with however c± = 1), and we have inserted a power of a, correspond- 
ing to the renormalisation constant Z discussed earlier, to make the result dimensionless. 

This shows that the (n — 1/n) dependence for the exponent of the Tsallis entropy is a 
general property of the continuum theory. Differentiating at n — 1, we find the main result 
of this section 

S A ~ -(c/6) log(mo) = (c/6) log(£/a) , (55) 

where £ is the correlation length. We re-emphasise that this result was obtained only for the 
scaling limit £ ^> a. However, for lattice integrable models, we shall show how it is possible 
to obtain the full dependence without this restriction. 

So far we have considered the simplest geometry in the which set A and its complement B 
are semi-infinite intervals. The more general case, when A is a union of disjoint intervals, is 
more difficult in the massive case. However it is still true that the entropy can be expressed 
in terms of the derivative at n — 1 of correlators of operators $ n . The above calculation can 
be thought of in terms of the one-point function ($ n ). In any quantum field theory a more 
general correlator (IliLi ®±n{ w i)) should obey cluster decomposition: that is, for separations 
\wi — Wj\ all 3> £, it should approach ($ n ) k . This suggests that, in this limit, the entropy 
should behave as Sa ~ A(c/6) log(£/a), where A = k is the number of boundary points 
between A and its complement. This would be the Id version of the area law(lQ(. When the 
interval lengths are of the order of £, we expect to see a complicated but universal scaling 
form for the cross-over. 
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B. Free bosonic field theory 



In this subsection we verify Eq. (155]) by an explicit calculation for a massive free field 
theory (Gaussian model.) The action 

S = j\ ((d^) 2 + m V) (56) 

is, as before, considered on a n-sheeted Riemann surface with one cut, which we arbitrarily 
fix on the real negative axis. 

To obtain the entanglement entropy we should know the ratio Z n /Z n , where Z n is the 
partition function in the n-sheeted geometry. There are several equivalent ways to calculate 
such partition function. In the following, we find easier to use the identity 3l|] 





^logZ n = -\f G n (r,r)<fr , (57) 

where G n (r, r') is the two-point correlation function in the n-sheeted geometry. Thus we 
need the combination G n — nG\. G n obeys 

(-V 2 + m 2 )G n (r,r') = 5 2 (r-r'). (58) 

Its solution (see the Appendix) may be expressed in polar coordinates as (here < r, r' < oo 
and < 9, 6' < 2im) 

4&./;« Ji/ -'S (V) tf/) , (59) 

where Ck{0,8') = cos(k6/n) cos(k6'/n) + sin(k6/n) sin(k8'/n), <i = 1, t4>o = 2, and Jk{x) 
are the Bessel functions of the first kind. At coincident points (i.e. r = r', 9 = 6'), and after 
integrating over 9 and A we have 

oo 

G n (r) = G n (r, r) = J2 d k I k / n {mr)K k/n (mr) , (60) 

fc=0 

where hix) and K k (x) are the modified Bessel functions of the first and second kind respec- 



tively [18]. 



The sum over k in (1601) is UV divergent. This reflects the usual short-distance divergence 
which would occur even in the plane. However, if we formally exchange the order of the sum 
and integration we find 
d 



dm 2 



I r \ roo | 

\ogZ n = - J d 2 rG n (r)=-^2d k J I k/n (mr)K k/n (mr)rdr = ^Yl d k k ■ ( 61 ) 

k " k 
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Interpreting the last sum as 2((— 1) = —1/6, we obtain the correct result, which we now 
derive more systematically. 

Let us first regularise each sum over k by inserting a function F(k/A n ): F is chosen so 
that -F(O) = 1, and all its derivatives at the origin vanish: however it goes to zero sufficiently 
fast at infinity. Since k/n is conjugate to the angle 9, we should think of this cut-off as being 
equivalent to a discretisation 59. Thus we should choose A n = A • n, where A ~ (59)~ l . 



To perform the sum over k, we use the Euler-MacLaurin (EML) sum formula 18] 

I 00 roo 1 00 R 

? E d„f{k) = f{k)dk - -f (0) - J2 T^/^CO) , (62) 



where I?2n are the Bernoulli numbers 18j. Using standard identities of the Bessel function 
(namely d k K k (x)\ k=0 = and d k I k (x)\ k=0 = -K (k) Q) ; we obtain 

G n (r) = 2j o dkI k/n {mr)K k/n {mr)F{k/nA) + —K%{mr)+J o rdr ^ —J^D 2j+1 {r) , 

(63) 

where we define D^x) = d l (I k (x)K k (x)) / dk l \ k= Q. In the last term in Eq. fl63l the order of 
the integral, derivative, and sum can be exchanged and each term in the sum is 

[)i roo Qi u 

— j o xdxI k (x)K k (x) = - — - = forz = 2i + l>2, (64) 



i.e. the sum in Eq. fl63l) is vanishing. 
Thus, still with the regulator in place 





dm 2 



roo roo | 

\ogZ n = -J^ rdr dkI k/n (mr)K k/n (mr)F(k/nA) - (65) 



where we have used J °° rK^imr^dr = l/(2m 2 ). 

Now the point is that in the integral over k the factor of n can be scaled out by letting 
k — ► nk. Thus this potentially divergent term cancels in the required combination G n — nG\. 
Having taken this combination, we may now remove the regulator to find the main result 



d \ogZ n 1 ( 1 
n 



(66) 



dm 2 Z n 24m 2 \ n, 

The integration of the last expression wrt m 2 (made following the recipe for the integration 
limits given in the previous section) gives 

Z n log a 2 m 2 ( 1 \ 
logTrp = log— = 24 [n- -) , (67) 
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and finally the entanglement entropy 



d 

S = — Trplogp = — —Tip" 

On 



d 



= (m 2 a 2 )*n n n, 

. On 

n=l 



12 

n=l 



-— logmV , (68) 



that agrees with the general formula we derive ( 1551) . with c = 1 and m = £ 
C. Integrable models and the corner transfer matrix 

In this subsection we verify ( |55l) for the transverse Ising chain and the uniaxial XXZ 
Heisenberg model. These results are also an independent check of the uniform convergence 
of the derivative wrt n of p n when n — > 1, since for these models we can compute the 
eigenvalues of pA exactly. 

Although these systems are integrable and their ground state is known, a direct calcula- 
tion of p is difficult. The difficulties arising in a direct calculation can be avoided mapping 
the quantum chains onto two-dimensional classical spin systems. As firstly pointed out by 
Nishino et al. [lij the density matrix of the quantum chain is the partition function of a 
two-dimensional strip with a cut perpendicular to it. In fact the ground state of a quantum 
chain described by an Hamiltonian H is also eigenstate of the transfer matrix T of a classical 
systems satisfying [H, T] = 0. Therefore the reduced density matrix of a subsystem A of 
the chain (defined in the Introduction as pa = Tr#|\I/) with B the complement of A) is 
the partition function of two half-infinite strips, one extending from — oo to and the other 
from +oo to 0, with the spins in B identified. 

This partition function is the product of four Baxter corner transfer matrices (CTMs)[2( 
A. If the lattice is choose in a clever way (i.e. rotated by 7r/4 with respect to the cut) one 
ends in the infinite length limit with [21] 

p A = A 4 = e - &c ™ , (69) 

where -£/ctm is an effective Hamiltonian, which, for the models under consideration, may be 



diagonalised by means of fermionisation (see for more details about this equivalence 21[ and 
references therein.) Note that TrpA 7^ 1, thus the usual density matrix is pa = Pa/T^Pa- 

The method just outlined is very general. However, for the integrable chains under 
consideration (and indeed for any model satisfying suitable Yang-Baxter equations [3]) it is 
possible to write H CTM = eO, with e the scale giving the distance between the energy levels, 
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and O is an operator with integer eigenvalues (for the Ising and XXZ models it is expressed 



in terms of free fermions.) Using this property the entropy is given by 32| 



e tv 1 t Pa log pA . . ^ . d\ogZ 

S = —Tip a log p A = -Tr ; h log Trp A = -e — h log Z , (70) 

TrpA 

where we defined Z = Tip a = Tre~ Hc ™. 

The Ising model in a transverse field can be described by the one-dimensional Hamiltonian 

tf/ = -£<-A£<7X + i, (7i) 

n=l n=l 

where a n are the Pauli matrices at the site n, and we normalise the Hamiltonian (following 



2 by imposing the transverse field in the x direction to be one, thus the transition is 
driven by the parameter A. The classical equivalent of flTTj) is two dimensional Ising model. 
For A = the ground state of the Hamiltonian (1711) is a quantum "paramagnet" with all 
the spins aligned with the magnetic field in the x direction, and (a^) = 0. In the opposite 
limit A = oo the magnetic field is negligible and the ground state is ferromagnetic with 
(cr*) = ±1. The (second-order) transition between this two regimes happens at A = 1. The 
exponent characterizing the divergence of the correlation length is v — 1, i.e. £ ~ |A — 
The CTM of the Ising model may be diagonalised in terms of fermionic operators. The 



CTM Hami 
and 1) is 



tonian, written in terms of the fermion occupation number hi (with eigenvalues 
2i] 



Hqtm = e A • ( 72 ) 

3=0 



The energy levels are 



( 2j + 1 e for A < 1 , K(y/T=W) 
€j = with e = Tr y \ 1 73 

I 2je forA>l, K(k) 

where K(k) is the complete elliptic integral of the first kind [lsj], and k = min[A, A -1 ]. 

For A < 1 

oo 

Z = Tre^ c ™ = [] [l + e- e(2j+1) ] , (74) 

3=0 

and the entropy from Eq. (I7D|) is 



oo 



s= ^ E T ^i^+I:'og(l+ e - ra+I, •). (75) 

3=0 1 + e j=0 



Analogously in the quantum ferromagnetic phase (A > 1) 



5-E7^ + EMl + e^). (76) 

3=0 1 ^ C j=0 
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FIG. 1: Entanglement entropy for the ID Ising chain as function of A. The dashed line is the limit 
for A — * oo, i.e. log 2. 

Fig. [1] shows a plot of the entropy as function of A, characterised by a divergence at the 
quantum critical point A = 1. Note that S(0) = and S(oo) = log 2, in agreement with the 
expectation that the pure ferromagnetic ground state (A = oo) has two possible accessible 
configurations with opposite sign of magnetisation (i.e. S(oo) = log 2) and the pure quantum 
paramagnetic ground state (A = 0) has only one configuration available with all the spins 
aligned in the direction of the magnetic field x and the resulting entropy is zero. 

Let us analyze in details the behaviour at the critical point. For A — > 1, e — ► in both 
the phases, thus the sums Eqs. (!75j) and ( J76l) can be approximated by the integral 

with ~ we mean in the critical region. The A dependence (we use K(0) = tt/2, K(x) = 
-1/2 log(l -x) + 0((1 - x)°) H) 

S ~ l2~7 " ~I2 l0g(1 ~ k) = h + Cl ' (78) 
where in the last equality £ oc |1 — k^ 1 has been used. This agrees with [55J, with c = 1/2 
for the Ising model. The constant C\ is not universal. 

Another model whose density matrix has been derived by using the CTM is the XXZ 
model 

H xxz = + « + i + A« + i) , (79) 

n 
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for A > 1 (whose classical equivalent is the Baxter six- vertex model [201] .) This model has 
an Ising-like ferromagnetic state for A > 1 and a planar XX ferromagnetic ground state at 
< A < 1 (the case A < describes the antiferromagnetic regime, in which we are not 
interested.) At A = 1 the Hamiltonian is isotropic (XXX) and approaching such a point 
from large A values the correlation length diverges as (see e.g. [20I ) 

because at A = 1 a massless excitation (Goldstone mode) is present in the spectrum. 
The CTM Hamiltonian Hqtm for this model has been obtained in Ref. 21 ] 

00 

Hqtm = £ 2 « > ( 81 ) 

with e = arccoshA, for A > 1. The entropy is given by Eq. (1761) . Close to the isotropic 
point (A = 1) it holds e ~ v^VA - 1 and (see Eq. (J77J)) 

(82) 

12 e 12v / 2 v / A^l 
Using Eq. (|80[) to write A in term of the correlation length, we have 

S^~log^, (83) 
o a 

in agreement with ESI with c = 1. Again, in the limit A — ► 00, S = log 2, since the 
ferromagnetic ground state is Ising-like. 

The method used here for the transverse Ising model and the XXZ model can in principle 
be applied to all those integrable model whose weights satisfy Yang-Baxter equations. 



V. FINITE-SIZE EFFECTS 



So far, we have studied the entropy either at the critical point, for a finite subsystem, 
or away from criticality in an infinite system. These two regimes may be linked by a 
generalisation of finite-size scaling theory. This would assert, for example, that the entropy 
of a subsystem A which forms e.g. the left half of a finite non-critical system of length 2L 
should have the form 

S4L,0 = ~(logL + s F ss(V0) , (84) 
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with sfss(O) = (this because we are referring to the scaling part, the constant term 
found before is not universal and it defines an overall additive normalisation of S) and 
s fss( x ) ~ — logx, for large x, so as to recover S = (c/6)log£ in the infinite length limit. 
sfss^) should admit a small x expansion 

sfss(^) = J2 s i x23 ' ( 85 ) 
with Sj universal coefficients, and a large x expansion 

s°° 

sfss{x) = -loga; + ^^-, (86) 
;>o xJ 

with s°° also universal (in special situations logarithmic corrections could also be generated.) 
In the next section we test this hypothesis for the free massive field theory. 



A. FSS in the Gaussian model 



Once again we consider the n-sheeted surface, but it now consists of n discs of finite 
radius L, sewn together along the negative real axis from — L to 0. The Gaussian two-point 
function in a finite geometry has been calculated in the Appendix Eq. ( 1A2I) . Setting 9 = 9', 
r = r' and integrating over 9' the propagator at coincident point is 

- °° 2/L 2 J 2 k/n (a k/n , t r/L) 
G n (r) = E^E 72 T— l . )a 2 JTjTT —2 » ( 87 ) 

where a V) i denotes the i— th zero of J u (x). G(r) is the analogous of (joOj) in a finite geometry. 

As we showed in the previous section, the right order to proceed to get a sensitive Z n is 
first to perform the sum over i, then sum over k and finally integrate over r. This is really 
hard, requiring a sum of Bessel functions for generic argument, over the zeroes of different 
Bessel functions. What one can do is inverting the order of the sums and integrations and 
try to understand what happens. From the previous exercise we know that this operation 
has to be done with care. 

However, the formal result is 

logZ„ = ^d tS log j^g±g . (88) 

Since the large % behaviour of the zeroes of the Bessel functions is a v ^ ~7r(i + z//2 — 1/4) 
(see e.g. [la]), this sum diverges. 
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A first, qualitatively correct expression for the universal function sfss(^) can be obtained 
from Eq. fl88l) . truncating the EML formula Eq. fl62l) in the variable k at the first order 
in the derivative. This approximation gives a correct form of the result because it takes 
into account the complete infinite volume result (whose truncated expression is exact.) The 
integral term in the EML is divergent, but it cancels as before in the ratio Z n /Z n if the 
cut-off in the angular modes is properly chosen as in the infinite volume case. The EML 
approximation at the first order is (we do not write the integral) 

_ 7rL 2 (q" 2 - m 2 ) ^ ^0(0^) _ F^mL, L/a) 

° g n ~ 12n ^(ooi + ^lWi + WW " 



where we used 

da k i 



dk 

which can be derived using the formulas 



2 71(00,,)' 1 J 

181 ] of the derivative of the Bessel functions with 



o 



respect to the order. The remaining sum over i (the zeroes of the Bessel function) is now 
finite. 

This gives the entropy as 



d 

S = - — Trf 

on 



2F 1 (mL,L/a), (91) 



n=l 



which reproduces the correct limits for m = and L = oo. 

From this formula we may compare with the FSS ansatz ( 1841) . with c = 1: 

s? ss (mL) = l2(F 1 (mL,L/a)-F 1 (0,L/a)) = -^ ^j JT+L^ ' < 92 ) 

that, as expected, is a function only of the product mL and all the dependence upon a 
disappears after the subtraction. This agrees with sfss(O) = and sfss(O) = — logx, for 
large x (it can be easily shown, since in this limit the sum can be replaced by an integral.) 
The over-script (1) is there to remind that we made a first-order approximation in EML. 

Eq. fl92|) characterises completely the crossover between the mass dominated (non- 
critical) and the geometry dominated (critical) regime. 

A plot of Spg S is shown in Fig. [2] (calculated as the sum of the first 1000 zeroes of Jo.) 
From this figure we see that an optimum approximation for x > 2 of Spg S (x) is Sfss( x ) = 
- \ogx+s^, with s§° = 0.120912 (this value of is a fit.) We can also calculate analytically 



all the universal coefficients of the small mL expansion 



22|: 



,f) = _,f^ = _l ( 93 ) 
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FIG. 2: Spgg(ic): Exact expression obtained as numerical sums over the first 1000 zeroes of the 
Bessel functions compared with small (quadratic, quartic and sextic) and large x (log cc + Sg°) 
approximations. Even the use of quadratic and large x approximation may reproduce the right 
formula over all the range. Inset: Comparison of several small x approximants for x < 2. 



b 2 



s 3 



b 4 



^Yo(a ,i) 1 



OO 



i Ji( a o,i) «o,-, 
Y (a 0t i) 1 
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J= i Ji{a ,i) al 
1.15 x 10" 3 



s 5 



3456 ' 

-1.98 x 10" 4 , 



s (i) 



3.42 x 10" 



(94) 

(95) 
(96) 



Fig. [2] (see also the inset) shows a plot of s 



(i) 

FSS 



X 



compared with quadratic, quartic, and 



sextic approximation. The agreement is excellent in the region x < 2. 

Before starting the full calculation of the function spss i x ) let us summarise what we can 
learn from the first order approximation in the EML expansion: A rather good approxima- 
tion of the full function may be obtained by matching only the first-order small argument 
expansion with the large x behaviour ~ log x (which is exact and does not depend upon the 
approximation.) Thus in the following we show how to calculate the coefficients Sj without 
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any approximation. 

Our starting point is again Eq. (1BB1 . but this time we will not use EML sum formula. 
For this reason we have to be care since the sum is divergent. 

The interesting object is \ogZ n — nlogZ. In particular we can write an FSS ansatz also 
for this "free-energy" and so subtract the m = part, obtaining the universal function 
f(x = mL) 



k i 



n log 1 + 



x 



at 



log 1 + 



X 



k,i / \ ®k/n,: 

. This formula has to be intended in a formal way: in fact 
it is the difference of two diverging sums. To make this difference sensible the recipe for the 
cut-off explained in the previous section has to be used. Anyway from the computational 
point of view it is simpler to make the calculation without care of the cut-off and to adjust 
only at the end the result. The numerical sum over the zeroes of the Bessel functions of 
generic order cannot be done as before, but the small x expansion of f(x) = fj x ^ with 
coefficients 



(97) 



form which Sfss(^) = — | x 



J k i 



n 



a 



2j 
k.i 



2j 
a k/n,i 



can be calculated analytically. For example 



wrong 



d_ 

dn 



n=l 



a j 



-1 + 



2k doth,- 
oi k ,i dk 



a 



k.i 



that using [22 1 



leads to 



wrong 







oo j 

i=i "fc,i 



k 



3 ^ dk dkA(k + l) 



4{k + l) 



1 



i(2^ 
4 V 6 



(98) 



(99) 



(100) 



101) 



We use the over-script "wrong" because it is not strictly correct to make the calculation this 
way, since we are implicitly using the same cut-off for Z and Z n . In fact, applying the EML 
to the sum in ( 110 II) . one has 



wrong 



3 V , 1 



4 I jo 



dk 1 

h -2 + 

{k + lf 6 



(102) 



that has a first order term = — 1/4 (in agreement with what previously found) but the 
integral is not vanishing. This is a finite difference of the two divergent expressions with the 
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wrong cut-off. The value of the integral (2 / dk(k + 1)~ 2 = 2) must be properly subtracted 
to have the right result 

3 7T 2 

Sl = — (2 3) = -0.217401... . (103) 

4 6 

This is close to the first order EML result —1/4, signaling that such approximation, not only 
reproduces the qualitative physics but is also quantitatively reliable (at the level of 10%.) 
In the same manner one can calculate all the coefficients Sj, using the more complicated 
expressions for the sum of higher negative powers of zeroes of Bessel functions reported in 
the literature 221. 



VI. HIGHER DIMENSIONS: SCALING OF ENTROPY AND AREA LAW 

The scaling hypothesis plays a fundamental rule in understanding classical phase tran- 
sitions. Crudely speaking, it asserts that the microscopic length scale a does not enter 
explicitly into thermodynamic relations near the critical point, as long as the various ther- 
modynamic variables are suitable normalised and allowed to scale with their non-trivial 
scaling dimensions, related to the various universal critical exponents. 

In fact, from the scaling of the singular part of the free energy density (here t = \T— T c \/T c , 
£ oc t~ u is the correlation length, h the external magnetic field, f±(x) a universal function, 
the subscript ± refers to the two phases, and is the scaling dimension of h) 

f Bi n S (t,h)=C d f±(hC yh ), (104) 

the critical behaviour of all the thermodynamic observables and in particular the scaling laws 
may be derived. Note that / s ; ng is not the total free energy density: there is another non- 
universal piece which has an explicit a~ d dependence on the microscopic cut-off. However, 
this term is analytic in the thermodynamic variables. 

As argued by Srednicki[10], for d > 2 the entropy Sa is proportional to the surface area 
A of the subsystem A. Thus we should discuss the entropy per unit area s = Sa/ A. In 
analogy with the classical case, we may conjecture a scaling form for the singular part of 
the entropy per unit area near a quantum phase transition 

s sing (g,h,T) = C {d - 1] s±(hC yh ,TC z ) , (105) 
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where yh > describes the relevant effect of the field conjugated to the order parameter, 
z > the relevance of the temperature close the the quantum phase transition (in all the 
examples considered up to now z = 1), and all other neglected operators are supposed to be 
irrelevant. s± should be a universal function apart from the normalisations of its arguments. 
The relation between £ and g (the parameter driving the transition) is £ = \g — g c \~ u - As for 
the free energy, s sing is not the total entropy: we also expect to find explicitly a-dependent 
pieces which are, however, analytic in g. 

From the scaling of the entropy all the scaling laws can be obtained, e.g. the specific heat 
goes like 

C = ~ T £-(d-i+z) _ ( 106 ) 

Similar identities can be derived for other observables. 

This scaling can be explicitly checked in the case of the Gaussian model in all dimensions 
for h = T = and z — 1. To calculate the entropy of a rf-dimensional system close to a 
quantum critical point, one has to consider a d + 1 field theory. We consider the geometry 
where d—1 dimensions are translational invariant in a domain delimited by an hypersurface 
of area A, and on the remaining two dimensional plane there is a branch cut going from 
to oo in an arbitrary direction. The choice of a different geometry is not expected to change 
the main results. 

The two-point function is translational invariant in d — 2 directions and coincides with 
Eq. (}59l) on the cut plane, i.e. 



G(rr, r, 0, r', V) = / g*^_L f * £° ^^W ^ ff) (m) 

where r±_ is the vector between the two considered points in the d—1 dimensional subspace. 

The calculation proceeds as in the one-dimensional case with the substitution m 2 — > 
m 2 + k\, and the integration over k± in front of all the equations. The space integration is 
obviously over d d ~ 1 r' 1 rdrd9 (we use r'j_ to stress that it is different from that appearing in 
Eq. (11071) that instead is r± = 0, since the two-point function at coincident point has to be 
considered.) Thus the analogue of Eq. (I6T1) is 

log Trp = / ^-Vl / p^.. \, (» - ") ■ (108) 

To have a finite result both the integrations are in a finite region. This means that the 
system must live in a finite box with / d d ~ l r'< L = A and a cut-off proportional to a -1 must 
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be understood for the integration over k±. 
In this way the entropy is 

A f d d -^k ± hj + m 2 
S = ~uJ (2^ l0g fciT^- (1 ° 9) 

In general, the integral diverges like a^^ 1 ^. However the coefficient of this divergence (and 
other sub-leading divergences which can occur for sufficiently large d) are in general analytic 
in m. There is, however, a finite piece which behaves as m d ~ 1 in agreement with (11051) . Note 
however that as d — > 1 this singular term combines with the non-singular piece oc a~^ d ~ x > to 
give the previously found result oc log(ma). 

Just as for the classical free energy, one may also conjecture a finite-size scaling form for 
the entropy of the form 



s 



sing 



[L, g, h, T) = L-«-Vs 1?ss {L 1/v \g - g c \,V h h, L Z T) . (110) 



For h — T — it reads 

s sing (L,g,h = T = 0) = A s L-( d -Vss(V0 , (HI) 

with sfss(£/0 a function which satisfies sfss(0) = 1 and sfss(^) oc 1/x for large x, in 
order to reproduce the infinite volume limit, and A s is a non universal constant fixing the 
normalisation of the entropy In one dimension, this relation again does not make sense. The 
logL behaviour may be seen as resulting from a cancellation between the L - ^ 1 ) behaviour 
of s a i n g against a constant term coming from the regular part, with A s = 0((d — 1) _1 ) as 
d — *■ 1. This gives the anticipated form 

S sing (L,£) = | (logL + WVO) , (H2) 
that has been explicitly checked for the Gaussian model in Sec. |V] 
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VII. NOTE ADDED. 



We recently were made aware, some time after this paper appeared in print, that the 
result (3.32) for more than one interval is incorrect in general. This error can be traced to 
the assumption at the beginning of Sec. 3C that there exists a uniformizing transformation 
for the n-sheeted Riemann surface into the sphere. In fact the surface has in general a 
non-trivial genus and therefore such a transformation does not exist. For the case of two 
intervals and n = 2 it can readily be seen that the surface is topologically a torus, and 
therefore the partition function Z 2 depends on the whole operator content of the theory, not 
only the central charge as claimed. For other values of n the partition function should also 
depend on the structure constants (OPE coefficients). It is therefore highly unlikely that 
this would not remain true for the derivative at n — 1. Indeed, a study of the computation 
for the 2-interval case for a compactified boson, first carried out in 1987 by Dixon et al. 24| , 
shows that this is the case, and that our (3.32) is incorrect. This is also the result of recent 
numerical investigations 25|, |2fjj]. We are grateful to V. Pasquier and F. Gliozzi for pointing 
out this work. 

Although this error also carries over to the case of a semi-infinite system, mentioned just 
below (3.32) and also to Eq. (3.33), it does not impact on our later application of these 
methods to studying the time-dependence of the entanglement entropy following a quantum 
quench 27| . This is because in this case only the limiting behavior, when the various cross- 



ratios are either small or large, is needed, and this is insensitive to the precise form of (3.32). 
Likewise, as far as we are aware, all the other conclusions of the present paper remain valid. 



APPENDIX A: TWO-POINT FUNCTION IN A n-SHEETED RIEMANN GE- 
OMETRY: INFINITE AND FINITE VOLUME RESULTS 

The Green function G(r, r') of the Helmholtz differential equation (— V 2 + m 2 )f(y) = 0, 
with specific boundary condition, is the solution of 

(-V^ + m 2 )G(r,r') = 5 d (r-r'). 

G(r, r') admits an eigenfunction expansion 

G(r,r / ) = E Ar ^( r )^( r ') ) 

k 
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in terms of 0&(r), eigenfunctions of the Helmholtz differential operator (— V 2 + m 2 ). N k is a 
normalisation constant that should be derived from the orthonormality requirement of the 
eigenfunctions 

N m / m (r)</> n (r')d^ = 5 nm . 



In the case of a n-sheeted Riemann surface we are interested in, the eigenvalue problem 
is solved in polar coordinates r = (x,y) = (r cos 9, r sin 9). A complete set of eigenfunctions 
is 

</>l ti = cosOVfl) J v (Xir) , (f> b U;i = sm(iv6)J v {\ir), (Al) 

with J u (x) Bessel functions of the first kind. Note that the Bessel functions of the second 
kind Y v (x), that are also eigenfunctions of the same differential operator, do not enter in the 
expansion since we require regularity at r = 0. 

Let us consider first the solution in a finite geometry. The infinite volume limit is recovered 
by taking the limit L — > oo in the sense of distributions. Imposing the 2irn periodicity 
boundary condition we have that v is an integer multiple of 1/n, i.e. v — k/n. Constraining 
the eigenfunctions to vanish at r = L, the eigenvalues are \{L = a u ^, with a v ^ the i-th zero 
of the Bessel function. 



Using the orthogonality relation of the Bessel functions (see e.g. |l8l |) 

L £2 

rdrJ v (a u>i r/L)J u '(a v ' ;i r/L) = — J v+X {a v>i )5 vv > , 

we get from the orthonormality requirement 

_ 4 



27m Jl /n+1 {a k/n ^ ' 
that gives the Green function 

G(r, e, r>, e<) = ^±d k ± J2 2 'f Jh/A^^ e >) . (A2 ) 

The factor d k comes from J 27r ldx = 2n and / 2 ^ dx cos 2 x = J 27r dx sin 2 x = 7r. 

In the limit L — > oo, the index z becomes continuous a u> i/L — > A and the in the 
orthonormality condition is replaced by 5(A — A'). The limit of the normalisation factor is 
-Wfc(A) = d k \/ (27m), that leads to the two-point function reported in the text (l59il . Note that 
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a slightly different form of this two-point function (satisfying Dirichlet boundary conditions) 
was used to investigate the critical behaviour at an edge [23]. 
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